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We study triplet pairing correlations induced in an SFS trilayer (where F is a ferromagnet and S 
an ordinary s-wave superconductor) by spin flip scattering at the interfaces. We derive and solve 
self consistently the appropriate Bogoliubov-de Gennes equations in the clean limit. We find that 
the spin flip scattering generates m = ±1 triplet correlations, odd in time. We study the general 
spatial behavior of these and of m — correlations as a function of position and of spin-flip strength, 
Hspin- We concentrate on the case where the ferromagnet is half-metallic. We find that for certain 
values of H sp i n , the triplet correlations pervade the magnetic layer and can penetrate deeply into the 
superconductor. The behavior we find depends very strongly on whether the singlet order parameter 
is in the or ir state, which must in turn be determined self-consistently. We also present results 
for the density of states (DOS) and for the local magnetization, which, due to spin-flip processes, 
is not in general aligned with the magnetization of the half metal, and near the interfaces, rotates 
as a function of position and H 3p i n . The average DOS in both F and S is shown to exhibit various 
subgap bound states positioned at energies that depend strongly on the particular junction state 
and the spin flip scattering strength. 

PACS numbers: 74.45.+C, 74.78.Fk, 74.78.Na 



I. INTRODUCTION 

Nanoscale structures involving Ferromagnet (F) and 
Superconductor (S) junctions illustrate the unique inter- 
action of superconducting and ferromagnetic symmetries 
and provides a novel opportunity to study the powerful 
influence that the spin degree of freedom plays in trans- 
port and thermodynamic properties of such systems. The 
now well established variety of phenomena 1 induced by 
the resulting proximity effects, includes exotic singlet su- 
perconducting correlations, in particular the damped os- 
cillatory Cooper pair amplitude in the magnet, with a 
spatial decay length in the clean limit over a few nanome- 
ters for strong magnets (such as Ni, Co, and Fe), and con- 
siderably less than the superconductor coherence length 
£o- These oscillations lead to the possibility of switch- 
ing between and ir junction states, with considerable 
potential 2 for applications. The superconductor region 
correspondingly becomes affected as it experiences the 
pair-breaking effects of the ferromagnet and becomes lo- 
cally magnetized. These mutual effects depend consid- 
erably on the strength of the magnet and transparency 
of the interfaces, usually assumed spin-independent. If 
the interface scattering is generalized to include spin de- 
pendence, where the spin of the impinging electron is 
flipped when traversing the corresponding interface, the 
whole picture can be modified, including the emergence 
or enhancement of exotic triplet states, that involve odd 
frequency or different time triplet correlations. 

Triplet pairing correlations can arise in ferromagnet 
and superconductor heterostructures involving supercon- 



ductors with a rotationally symmetric pairing symmetry 
(s-wave pairing), since they involve odd time symmetry 
pairing, as originally proposed 3 in a different context. 
In F-S structures the proximity effects associated with 
the magnet break spin rotation invariance. The super- 
conducting order parameter in this scenario changes sign 
under time coordinate interchange of the two electrons 
comprising Cooper pairs. When a single quantization 
axis exists for the system, the only possible triplet pair- 
ing state is the one comprised of opposite spin pairs (the 
772 — projection on the given quantization axis). If there 
exists non-collinear or inhomogeneous magnetization in 
the system, as can occur in structures involving differ- 
ently oriented F layers or an in-plane spin flip scattering 
potential, equal spin m = ±1 triplet correlations can also 
arise. Several investigations into triplet effects in super- 
conductor and ferromagnet hybrids has revealed a host 
of interesting and exotic phenomena^^ including the 
possibility of a long-ranged superconductivity proximity 
effect in F-S structures. A new superconducting state 
that can potentially extend superconducting correlations 
into the magnetic region over long distances brings with it 
a host of useful device application involving low temper- 
ature nanodevices, including nanoelectromechanical sys- 
tems (NEMS), and superconducting circuits (ir junctions 
with J < 0). 

While there has been recently considerable interest in 
trying to isolate and detect the triplet pairing state that 
is predicted to exist in such S-F structures, it can be dif- 
ficult to disentangle the triplet and singlet correlations. 
It is therefore of interest to investigate heterostructures 
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that restrict the singlet order parameter somewhat, yet 
retain the desired triplet correlations. The pinpointing 
of triplet effects can be exploited with the use of highly 
polarized materials, namely half metallic ferromagnets, 
where only a single spin channel is present at the Fermi 
level. The ordinary singlet pair amplitude is thus sup- 
pressed, since the magnet behaves essentially as an in- 
sulator for the opposite spin band. Several half metallic 
candidates are considered in connection with supercon- 
ducting hybrids and spintronic applications. These in- 
clude the conducting ferromagnet£ CrC>2, the manganese 
perovskite££ La2/3Cai/3MnC>3, and the Heusler alloys, 9 
possibly Co2FeSi and CoaMnSi, which are attractive from 
a nanofabrication standpoint, since growth by sputtering 
techniques is applicable. This is pertinent since CrC>2 
cannot be grown by sputtering and is metastable. Thus it 
is of interest to determine the circumstances under which 
triplet correlations and related single particle signatures 
emerge when a wide variety of spin flipping strengths for 
a half metallic ferromagnet in F-S nanojunctions, with S 
being a conventional s-wave superconductor and interfa- 
cial spin flip scattering providing the required symmetry 
breaking. 

The spin flip processes and their participation in prox- 
imity effects have been explored in several different exper- 
imental setups 6 - ' 10 ' 11 ' 12 possible in many instances due to 
advanced e-beam lithography and sputtering techniques. 
There is no current experiment however, that offers com- 
pletely indisputable evidence of triplet correlations in fer- 
romagnet and conventional superconductor hybrids, and 
thus further experimental investigations are needed. The 
spin flipping associated with the intrinsic exchange field 
in the ferromagnet in S-F bilayers was linked to critical 
temperature variations^ - as a function of F layer thick- 
nesses. It was suggested that the measured Josephson 
current 6 in a sample with two NbTiN (s-wave) super- 
conductors coupled by half metallic CrC>2 is due to a 
supercurrent carried by spin triplet pairs since the elec- 
tronic transport in CrC>2 is metallic solely for the spin 
up band, and the expected magnet thickness for that 
system exceeded the estimated singlet correlation length 
for that sample. Current- voltage measurements jii in sin- 
glet superconductor-half metallic point contacts revealed 
a marked resistance decrease and increased normal state 
conductance at small voltages, attributed to spin singlet- 
triplet conversion. Single particle spectroscopy results for 
the density of states (DOS) were also reported^ 2 - for F-S 
bilayers including the strong ferromagnet, Ni. In that ex- 
periment, the conductance signature was measured as a 
function of ferromagnet thickness, dp, revealing an inter- 
esting double peak structure and other subgap features 
that could not be theoretically accounted for within the 
dirty limit framework. 

Numerous theoretical approaches involving spin de- 
pendent scattering of some sort have helped pave the 
way towards unveiling the role of triplet pairing corre- 
lations in diffusive SFS hybrid nanostructures or clean 
SFS junctions, both within the quasiclassical=. regime. 



The intermediate regime, separating diffusive and ballis- 
tic motion was also studied quasiclassicallyJ^ Any purely 
microscopic approaches, that retain quasiparticle infor- 
mation at the atomic scale, typically involved spin in- 
dependent scattering potentials at the interfaces ! 15 ' 16 If 
the pair-breaking mechanism of spin-flip scattering at the 
interfaces is included, the resultant interchange of spins 
yields complicated normal and Andreev reflection events. 
The investigation into these issues have been predomi- 
nately in the diffusive regime however. For instance, with 
relatively thick half metals, Josephson coupling can oc- 
cur via triplet correlations from the singlet superconduc- 
tor and spin mixing occurring at a spin active interface. 13 
For calculations neglecting the mutual influence of super- 
conducting and ferromagnet order parameters, spin flip 
scattering was shown to also have a detrimental effect on 
the residual supercurrent, thus limiting such junctions 
as useful spin switches^ A long range triplet compo- 
nent can arisei^ when the ferromagnet has a Neel do- 
main structure where the in-plane magnetization rotates 
with changing depth in the magnet. The odd frequency 
pairs arising from spin flipping at the junction interfaces 
can cause a peak in the local DOS of a diffusive half 
metallic ferromagnet^ If the ferromagnet can be mod- 
eled by a conical magnetic structure^ as in Holmium, it 
was found that both singlet and triplet correlations un- 
dergo short range decay. Also, the decay length of the 
Josephson current was shown to decrease 2 ^ with spin- 
flip and spin-orbit scattering, with spin-orbit scattering 
typically being the more destructive of the two. By il- 
luminating the junction with microwave radiation at the 
proper resonance frequency however, the critical current 
can be enhanced^ 2 - due in part to singlet-triplet conver- 
sion processes. Nearly all of the cases studied thus far 
involve the quasiclassical method, and it is unclear how 
this landscape is modified when atomic length scales are 
not eliminated in the pertinent equations and when self 
consistency of the singlet order parameter is taken into 
account. 

In this paper, we address some of the above issues by 
presenting a fully self consistent framework for a clean 
nanoscale trilayer junction comprised of a half metal 
sandwiched between two conventional, s-wave supercon- 
ductors. The pair breaking mechanism is spin flip scat- 
tering at the interfaces, which produces m = ±1 odd 
time pairs, and modifies the triplet m = component. 
The presence of two S layers, coupled through F via the 
proximity effect allows us to compare and contrast the 
and 7r states. Our method is based on the quantum me- 
chanical Bogoliubov-de Gennes (BdG) equations in the 
clean limit, which is ideal for half metallic ferromagnets 
and proximity effects that can involve singlet correlations 
in the magnet with very short decay lengths of just a 
few nanometers. We are able to fully take into account 
proximity effects in the magnet and "inverse proximity 
effects" that arise in the superconductor regions, includ- 
ing the presence of a magnetic moment component nor- 
mal to the magnetization in F. We employ a recently 
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developed 23 method to determine the triplet correlations 
in such structures using a Heisenberg representation to 
derive the time dependent quasiparticle wavefunctions. 
We consider spin-active interfaces by incorporating in- 
plane spin dependent scattering in the effective Hamilto- 
nian. By varying the spin scattering strength parameter, 
H S pim over a broad range, long range triplet correlations 
are shown to emerge and evolve. We study the spatial 
profile of all possible triplet correlations which depends 
on their corresponding projection onto the axis of quan- 
tization, which is taken to be along the fixed direction of 
magnetization in the half metal (the z axis in our case) . 
The relative admixtures of triplet amplitudes with to- 
tal spin projection m = on the z-axis (labeled /o) and 
those with m = 1 total spin projection quantum number, 
/i, depend crucially on whether the junction is in a or n 
state. The junction state with the lowest free energy, and 
its corresponding stability, is dictated not only by the ge- 
ometry and spin splitting strength of the magnet^ but 
also by the magnitude of spin scattering at the interfaces. 
An accurate determination of this requires a self consis- 
tent calculation: this is even more evident when one con- 
siders that the triplet correlations in general peak near 
the interfaces, where self-consistency is most critical. Af- 
ter presenting the triplet correlations within the system 
as a function of the spin scattering strength, we turn our 
attention to the effect spin active interfaces have on other 
physically important single-particle quantities such as the 
average DOS, and local magnetic moment. We find that 
the DOS has a subgap signature that depends on whether 
the junction is in a or tt state and the degree of spin 
scattering at the interfaces. 

The rest of this paper is organized as follows: in Sec. Hi] 
we discuss the methods we use to evaluate the order pa- 
rameter in a fully self consistent way and to determine the 
triplet correlations. As pointed out in Ref . fill self consis- 
tency is absolutely essential in order to correctly obtain 
the odd parity triplet correlations: without self consis- 
tency the Pauli principle is violated. We review and dis- 
cuss the appropriate quasiparticle expansions, the evalu- 
ation of the matrix elements, and other relevant details 



The first term in brackets is the single particle Hamil- 
tonian for a quasiparticle with effective mass, m, Fermi 
energy, Ep, and scattering from a spin independent po- 
tential Voir). The pair potential, A(r), characterizes the 
spatial dependence to the superconducting singlet corre- 



in the solution of the corresponding eigenvalue problem. 
The definition of the time-dependent triplet amplitudes 
is given and other quantities that are also of interest, 
such as the local magnetic moment and the local DOS 
are defined. Then, in Sec. IIII1 the condensation energy, 
which reveals the relative stability of the and tt phases, 
is calculated as a function of H sp i n . We then discuss in 
detail our results for the spatial and time behavior of the 
triplet amplitudes as a function of the spin flip scatter- 
ing strength. The associated penetration depths of the 
equal spin triplet amplitudes into the superconductor re- 
veal the long range nature of these correlations. Next, we 
focus on the averaged DOS in each region of competing 
order parameter symmetries, and display the spatial de- 
pendence of the local magnetic moment, as it relates to 
the inverse proximity effect. We show that the induced 
magnetic moment vector rotates near the interfaces. Fi- 
nally, we give a summary of our results in Sec. IIVI 



II. METHODS 

The SFS junction that we study is a trilayer structure 
infinite in the plane parallel to the interfaces, which we 
label the x — z plane, and with total length d in the y 
direction, normal to the interfaces. The width of each 
of the two superconductor layers is labeled by ds and 
that of the ferromagnet by dp. The superconductors are 
s-wave and identical. The entire structure occupies the 
space < y < d, with one superconductor occupying, 
< y < ds , the ferromagnet: ds < y < ds + dp, and the 
other superconductor, ds + dp < y < d. There is spin- 
flip scattering at the interfaces, which will be described 
below. 

To determine the equations governing triplet effects in 
SFS nanojunctions, we start with a fundamental quan- 
tity, the effective Hamiltonian, Ti. e s, written in terms of 
creation and annihilation field operators and vector Pauli 
spin matrices, er, 



(2.1) 



lations, and will be calculated in a self consistent fashion 
as described below. The ferromagnetic exchange field, 
h(y) = hoz, representing the ferromagnetism, is taken as 
constant in the F layer and vanishing in the two S layers, 
and it is along the z axis of quantization. This intrin- 
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sic exchange field in the magnet, favoring a given spin, 
thus contributes to the overall behavior of triplet correla- 
tions. The important spin flip scattering will be assumed 
to be confined to the two interfaces near y — d$ and 
y = ds + dp. It takes place in the invariant x — z plane: 
V • er = V x (y)a x + V z (y)a z . Its z component represents 
a less important local modification of the ho field, while 
V x is the spin flip term. We have taken V y = because 
of the geometry and also for convenience: including a v 
terms precludes the use of exclusively real numbers in the 
numerical diagonalizations and leads to additional tech- 
nical irrelevant complications. Each of the triplet states 
can potentially exist over large length scales, thus allow- 
ing competing orderings to coexist. 

To solve the problem we diagonalize H c s via a Bogoli- 
ubov transformation. The details are given elsewhere^ 
and need not be repeated here. Through the use of 
standard commutation relations, we end up after some 
straightforward algebra, a general coupled four compo- 
nent set of equations. This leads to a generalization of 
the textbook^ Bogolioubov-de Gennes (BdG) equations, 
which give rise ultimately to spin singlet and triplet am- 
plitudes. By making use of the Pauli spin matrices and 
of a set of Pauli-like matrices p in particle-hole space, the 



general time and spin-dependent BdG equations can be 
expressed compactly as, 



[ Pz ® (H i - [h z - V z )a z ) + (A(y)p x + V x l) ® a x 



where the four component wavefunction, 
a vector of quasiparticle amplitudes, $ 



]$n(y,t) 

(2.2) 



.9$, 



is 



(y,t) = 

(u n ](y),u nl {y),v n - ] (y),v nl (y)) T e le "*, where the su- 
perindex denotes transposing, the u na and v na have their 
standard^ meaning as quasiparticle amplitudes and e„ 
is the eigenenergy. We have assumed here that Ep is 
the same throughout the sample: the majority (+) and 
minority (— ) bandwidths in F are Ep ± ho- The single 
particle quasi one-dimensional Hamiltonian Ho becomes 
in our geometry 



Ho^— — 

2m dy 2 



■e x -E F + V Q (y), 



(2.3) 



where e± is the energy in the transverse direction. Car- 
rying the time derivative through, and taking the outer 
product in Eq. (|2.2p , we can rewrite Eq. (|2.2p in the much 
less compact but intuitively more immediate form: 
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(2.4) 



where the spin dependent interface scattering potential 
should be understood to be given in terms of delta func- 
tion scatterers: Vi(y) — Vi[5(y — ds) + S(y — (ds + dp))], 
and as explained above, i = x,z. The convenient di- 
mensionless parameter H spin = 2mV x /kp characterizes 
the strength of the interface scattering. The spin flip 
cc-component, V x (y), technically complicates the calcula- 
tion and prevents the simple splitting of the BdG equa- 
tions into two separate equations by means of symme- 
try relations, as in the case of collincar magnetizations, 
or when a single quantization axis exists for the whole 
system. Here, all four components are needed since the 
exchange field in the ferromagnet as well as the spin-flip 
potential break the spin rotation invariance. 



tential, valid for all temperatures, T, is given by, 



A(y) = 9 -^fT,H(y) v n(y) +"i(v»ttv)] tanh (§0' 

n 

(2.5) 

where the sum is over eigenstates (the index n now sub- 
sumes not only the quantized index in Eq. (|2.4j) but also 
the transverse energies ej_) which is performed over all 
eigenstates with positive energies smaller than or equal 
to the "Debye" characteristic energy cutoff u>d, and g(y) 
is the superconducting coupling parameter that is a con- 
stant go in the intrinsically superconducting regions and 
zero elsewhere. 



The triplet correlation functions, odd in time, which 
The general expression for the self consistent pair po- are the main subject of our study are defined 16 in terms 
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of the usual field operators as, 

/o(r,i) ee i[M(r,t)^(r,0)) + (^ 4 (r,i)^(r,0)>] 

(2.6a) 

/i(r,t) ee i[(^(r,t)^(r,0)> - <^(r,t)^(r,0)>], 

(2.6b) 

where we are clearly free to choose one time coordinate 
to be zero, without loss of generality. These correlation 
functions must be odd in time because of the Pauli prin- 
ciple. Hence they vanish identically at t = 0. At T = 0, 
these expressions are conveniently written in terms of the 
quasiparticle amplitudes 2 ^: 

fa(y,t) = ^^2[uni(y)vni(y) -u n i(y)v n ^(y)]e- Knt , 

n 

(2.7a) 

h{y,t) = ^^2[u n r(y) v nt(y) + u ni(y)vni(y)]e~ ient , 

n 

(2.7b) 

where all positive energy states are in general summed 
over. In practice, we find that at finite times, results be- 
come cutoff independent beyond a value a few ui p, . How- 
ever, to ensure the vanishing of the triplet components 
at t = 0, it is necessary to sum over a much larger energy 
range. 

Besides the pair potential and the triplet amplitudes, 
we can also determine various physically relevant single- 
particle quantities. One such important quantity is the 
local magnetization, which is a measure of the so-called 
inverse proximity effect, and can be particularly useful 
in characterizing the magnetizing effects in the super- 
conductor as a result of the localized spin flip interface 
scattering and intrinsic exchange field of the magnet. 
It can also serve as an effective self-consistent measure 
of the magnetization field in the half metal. Recent 
magneto-optical Kerr effect measurements 2 ^ of a super- 
conductor/ferromagnet bilayer revealed that the super- 
conductor became magnetized, illustrating the need to 
determine the spatial behavior of the magnetization fully. 
In the presence of spin flip scattering, the local magnetic 
moment m will depend on the coordinate y (in our geom- 
etry) and, in the presence of the spin flip term it will have 



in general both x and z components, m = (m x ,0, Tn z ). 
In terms of the quasiparticle amplitudes calculated from 
the self-consistent BdG equations we have, 

m z {y) = -n B J2[v^{y) - v 2 ni {y)], (2.8a) 

n 

m x (y) = -2fj, B ^2v n1 (y)v nl (y), (2.8b) 

n 

where n B is the Bohr magneton. The sums in Eqs. (|2.8p 
involve a sum over eigenstates, as in Eqs. (|2.5p and (|2.7[) . 
although the energies e„ do not now appear explicitly. 

A very useful tool in the study of these phenomena is 
experimental tunneling experiments, where spectroscopic 
information, measured ideally by an STM, can reveal the 
local DOS. Therefore we have computed here also the 
local DOS N(y, e) as a function of y. We have N(y, e) ee 
Ni(y,e) + N l (y,e), where, 

N ( r(y,e) = J2[uL(y)S(e-e n ) + vl <T (y)d(e + e n )}, a =].,[■ 

n 

(2.9) 

In order to numerically solve the problem, we must re- 
express the equations in terms of matrix elements in an 
appropriate basis. These matrix elements are obtained 
via projection upon a orthonormal complete set, that in- 
herently satisfies the boundary conditions of vanishing 
wavefunction at the outer edges of the trilayer structure. 
Thus we write u na (y) = \/tydY,q=iU™ q sm(qTry/d), 

Vna(y) = \/^X^i<g sm (9 7r y/ c 0> witn CT =T, I- in- 
serting these into Eq. (|2.4[) . we have the general AN x AN 
matrix consisting of a 4 x 4 array of block submatrices, 
each of rank N: 
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v x 


H~ 


V 
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-H+ 


v r 


\v 





V 


-n 



*„ = e„*„, (2.10) 



where we measure all energies in terms of Ep, 
so that e„ = e n /E F , Z Bz = mV z /(k F d), 
and the vector ty n is the transpose of 

( M nl' ' U nN> U nV ' ' ' ) U nN' V nli ' V nN> V nl ' " J "nw)' 

We find, after lengthy but elementary algebra the matrix 
elements: 



( i-k \ 2 e± 

fe) + E~ F - lTl 



d + 21 



K, 



±Z Bz (U^U^ 



]±Z Bz {TJ?>US 
I 



L J 



(2.11) 



where / = h^jEp and Z Bz = 2V Z /Epd. The important spin flip component, off the main diagonal, giving rise to 
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equal spin triplet correlations is, 

V* = Z Bw [uPuf ] + Ul 2] uf\ (2.12) 

where — sm{q-Kds /d), — sin(qTr(dp + ds)/d), 

K q 1] = uP/iqn), K [ q 2) = U^ 2) /(qtt), and Z Bx which 
arises from the spin-flip scattering is Z Bx = 2V X /Epd. 
This is related to the parameter H sp i„ = 2mV x /kp 
defined earlier to characterize the spin flip scattering, 
by H sp in = kpdZsx/2. Without loss of generality 
we can take V z = V x = V sp i n . The matrix ele- 
ments corresponding to the pair potential are, T>ij = 

2/{Epd) Jq dyA(y) sui(iiry/d) sm(jiry /d), recalling that 
A(y) vanishes in the magnet layer, due to the coupling 
g{y). Since we are not permitted to use previous sym- 
metry relations among the quasiparticle amplitudes and 
energies that reduced the matrix eigensystem to 2N, we 
are forced, as mentioned above, to solving the AN x AN 
system, and retaining only the positive energy states. As 
in previous work, the diagonalization is performed iter- 
atively until the self-consistency condition Eq. (|2.5p is 
satisfied. 



III. RESULTS 

The results of our calculations are described in detail 
in this section. We will measure all the lengths in units 
of the Fermi wave vector kp, and define the relative di- 
mensionless coordinate Y = kp(y — d/2), i.e. Y = is 
at the center of the junction. All times will be given in 
units of ujp 1 via the dimensionless time r = toot- In 
the ferromagnet we have for the spin up and spin down 
band widths, = Ep + h and E^ = Ep — h. The 
dimensionless measure of the intrinsic exchange energy 
in the magnet is / = h/Ep. All results below are for 
the half- metallic limit, 1=1, discussed recently in the 
context of spintronic materials.— The spin flip scatter- 
ing at the interface is characterized, as previously ex- 
plained, by the dimensionless parameter H sp i n for which 
we will consider values between zero and unity. Geomet- 
rically, we will consider a system consisting of two thick 
superconducting layers, each of a thickness ds such that 
Ds = kpds = 300. Since odd-frequency or different-time 
triplet states arise from magnetic effects at the interfaces 
and exchange field in the half metal, this broad range of 
parameters will give revealing hints as to their existence. 
The chosen ds considerably exceeds the superconduct- 
ing coherence length £o, which we take to be £o = 50k p 1 . 
Thus ds = 6£o , so that we can disentangle quantum inter- 
ference effects from the modified Andrccv and scattering 
events at a spin flip interface. The two S layers are sep- 
arated by a ferromagnetic layer, which must be taken to 
be thin enough so that the two superconductors are still 
coupled through the F material via the proximity effect. 
We take Dp = kpdp = 10. All results are computed in 
the low temperature limit. 
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FIG. 1: (Color online) The dimensionless condensation en- 
ergy (free energy at T = in units of A r (0)Ao, see text), 
versus the dimensionless parameter H sp i n characterizing the 
spin flip strength at the interface, for an SFS junction with 
a half metallic ferromagnet. For the geometry chosen, we see 
that the both the n and the states are stable for all val- 
ues of H sp i n considered and that the n state has the lower 
condensation energy except at small values of H sp i n . 



Each junction between two consecutive S layers can be 
of the "0" type (with the order parameter in both S layers 
having the same sign) or of the 'V type (opposite sign). 
The characteristics of a or tt junction are directly con- 
nected to the spatial behavior of the pair potential A(y) 
and, to determine its precise form, this quantity must be 
calculated self-consistently so that the resulting singlet 
pair amplitude corresponds to a minimum in the free en- 
ergy. The relative stability of the different states that 
may be obtained through self-consistent solution of the 
BdG equations depends on the free energy of the junc- 
tions. We therefore consider first the stability of the sys- 
tem for our parameter values and geometry. In Fig. [T] we 
plot the condensation energy (the free energy at T = 0) 
of the system in dimensionless form, that is, in units of 
A(0)Ao where N(0) is the usual single spin density of 
states in the S material, and Ao the bulk value of the gap. 
Thus, the quantity plotted would be —1/2 for a bulk S 
sample. This energy is calculated with high precision as 
explained in previous^ work. We see in the figure that, 
consistent with the results of Ref. [HI, the condensation 
energies are reduced, in all cases, from what they would 
be for a bulk S sample. For all values of H sp i n both the 
and tt configurations of the structure are at least locally 
stable, but in general non-degenerate, showing that in- 
deed the two S slabs are indeed coupled via the proximity 
effect. At very small values of the spin flip parameter, 
the configuration is the stable one, but this changes as 
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H sp in increases: there is a first order phase transition at 
H sp i n < 0.2 and at larger values of H spin the it configu- 
ration is the stable one and the configuration is much 
less stable. The configuration metastable minimum is 
shallowest at H sp i n « 0.8 where the condensation free 
energy has a sharp maximum. The condensation free en- 
ergy of the 7r state is more weakly dependent on H sp i n 
with a maximum near H sp i n = 0.82 much shallower than 
that found in the state. Note that by decreasing the 
width of S, the qualitative results remain, but the overall 
results are shifted towards zero, resulting in the possi- 
ble elimination of the zero state altogether. In general, 
computational convergence time is increased as the con- 
densation for a given state approaches zero. 

We turn next to the spatial dependence of the general 
complex triplet pairing functions f (Y, r) and fi (Y, r) as 
defined in Eqs. (|2.7[) . In Fig. Owe plot the corresponding- 
triplet amplitudes for each of the two types of solution 
(0 or 7r) as a function of the dimensionless coordinate Y 
and at fixed time, r = 20, for 6 equally spaced values 
of H SV in in the range < H spin < 1 (see legends). All 
amplitudes plotted are normalized to Ao/g, so that, if 
the similarly normalized ordinary singlet amplitude were 
plotted, it would reach unity deep in bulk S material. 
The value r = 20 is chosen as being near— that which 
maximizes the correlations, and is such that the triplet 
pairing states have penetrated most of the two supercon- 
ductor regions. Results for the first group of 4 plots are 
for the 0-state solutions, and with the real and imaginary 
triplet amplitudes labeled accordingly. The bottom series 
of four panels are for the ir junction counterparts. The 
range of Y included in the plots is, for clarity, somewhat 
narrower than the sample size: regions where the ampli- 
tudes are very small or zero arc omitted. One can see 
that the amplitude fi vanishes identically in the absence 
of spin-flip scattering, since in that case both the total 
spin and its z component are good quantum numbers. 
For finite values of the spin flip parameter, all possible 
projections of the total spin exist. The spatial symmetry 
of the singlet Cooper pair is also reflected in the triplet 
pairing states: it is evident from this and the next figure 
that if the singlet order parameter is in a junction state, 
the corresponding triplet amplitudes maintain that sym- 
metry. This holds true for the spatially antisymmetric 7r 
junction results as well (bottom 4 panels). 

Turning our attention to the real part of fx for the 
0-state we see that it shows a monotonic decline in mag- 
nitude from the interface, over about three to four co- 
herence lengths, for the largest two spin flip strengths 
(with a superposition of rapid oscillations) . The remain- 
ing weaker scattering strengths are quite different in that 
they yield nonmonotonic behavior with a maximum deep 
into S at about 2£o , and then decaying to zero at roughly 
4£ (hence for H spin = 0.8, this correlation dies about £o 
earlier). These amplitudes , Re/i, are predominately pos- 
itive for higher spin transparency (smaller H spin ) junc- 
tions, and then undergo a sign flip for the stronger H sp i n . 
If we examine now the imaginary component Im/i, still 



for the state, we see similar opposite parity effects sep- 
arating the strongest H sp i n from the weaker values, here 
however there is a clearer separation between the curves. 
The time dependence here is noticeably different than for 
the real part; the triplet correlations have a faster rate of 
propagation, in that they have reached deeper within the 
sample for the same r. We have emphasized the triplet 
amplitudes in the S region, however these plots reveal 
that besides the expected fact that (ip-\(y,t)tp-\(y,0)) is 
not destroyed by the half metal, by including the proxim- 
ity effects in a self consistent way, we found non- negligible 
different spin triplet pairing in the half metal, but with 
a smaller magnitude that f\. 

The 7r state results for the equal spin pairing correla- 
tions have markedly different profiles than those for the 
state: Besides being highly peaked at the interface, where 
spin-flip scattering originates, Re/i has a very weak de- 
pendence on H sp i n , with an abrupt emergence only for 
the highest H sp i n and then still decaying over roughly 
the same distance in S. Overall, the it state fi amplitudes 
are suppressed, even for the case when both the and 7r 
states have the same condensation energy [H sp in ~ 0.2). 
The diminished 7r state results arises partly from the sym- 
metry requirements imposed upon ff. fi(—Y) = — /»(Y), 
thus in F, the fi amplitudes vanish at Y = 0, which can 
constrain the overall longer range spatial behavior. It can 
be concluded that the singlet Cooper pair order param- 
eter that minimizes the free energy, and which |A(Y) is 
typically larger, does not necessarily result in the larger 
triplet amplitudes. The imaginary parts can be discerned 
from the figure, but clearly the imaginary parts to the 
amplitudes also do not simply differ by a phase. Also, 
some amplitudes oscillate (with periods on the order of 
£o) while other plots show simple declines, with atomic 
scale oscillations superimposed on the general profile. 

The m = triplet amplitude, /o, with zero projec- 
tion of the z-component of total spin, does not vanish, at 
any finite time, even for H sp i n = since the total spin 
(as opposed to its z component) is not a good quantum 
number in the presence of the F material. It does still 
vanish at t = because of the Pauli principle. Thus, 
in general, when a single quantization axis exists for the 
system, /o coexists with the ordinary singlet s-wave com- 
ponent. We expect the results for f to be different those 
for fi since f does not emerge solely from H sp i n '- there 
are two competing spin flip effects in the z direction: the 
magnetization of the half metal, and the spin dependent 
scattering at the interface. It is clear from the form of 
the BdG equations (|2.4p that these two effects compete 
against each other. For both the and 7r state, the abso- 
lute value of a given component of /o at fixed time in S is 
again a non-monotonic function of H sp i n , but the over- 
all dependence is visibly different. The maximum value 
of fo in S is always at the interface: this is as expected 
from the above considerations regarding quantum num- 
bers. The behavior with H sp i n echoes (at larger values) 
that found for f±, and for the same reasons. In all cases, 
these triplet correlations clearly pervade the thin F layer, 
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FIG. 2: (Color online) The /i and /o triplet pair amplitudes (Eq. (12. 7[) ) for a 0-junction (top 4 panels) and 7r-junction (bottom 
4 panels) plotted as a function of the dimensionless coordinate Y for several values of H ap i n as indicated in the legend. The left 
panels show the real parts while the right ones show the imaginary parts, for values of H sp i n ranging from to 1. All results 
are at a fixed value of the dimensionless time, r = 20. 



while their penetration into S increases only weakly with 

In Fig. [3] we show results for the triplet amplitudes 
in the same format as in Fig. [5] but at fixed H sp i n and 
several values of r, so that the explicit time evolution 



of the different-time triplet states can be visualized. As 
before, the real or imaginary part of fo or /i are appro- 
priately labeled whether discussing the or tt state con- 
figurations (top four or bottom four panels respectively) . 
All results are for an intermediate spin flip transparency 
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FIG. 3: (Color online) The temporal dependence of the and n state junction triplet amplitudes. They are arranged as in 
Fig. [2] and plotted as a function of Y for several equally spaced values of r. These results are at a fixed value of the spin flip 
parameter, H sp i n = 0.2, where both the -k state and the state are stable (see Fig. [TJ. 



H sp in = 0.2, where both junction states have very ap- 
proximately the same condensation energy (see Fig. [l}. 
The region of the sample shown here is wider than that 
in Fig. [5] because the spatial range over which the corre- 
lations extend is now wider, at larger times. We see than 
the triplet correlations, which of course vanish identically 



at r = already pervade the S layer at the earliest times 
shown. 

The real part of the amplitude /i has, in the 0-state, 
a maximum in the S region that keeps propagating out- 
wards as r increases, reflecting the longer penetration 
of the correlations into S. This maximum becomes shal- 
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lower with increasing r however. This increased pene- 
tration occurs also for both the real and imaginary com- 
ponents of /o, although in the latter case the maximum 
value of the amplitude in S occurs near the interface ex- 
cept in some instances at the longest times studied. The 
largest r studied was determined by the need to avoid fi- 
nite size effects: after the different-time correlations have 
pervaded the entire S portion of the sample the results 
would be contaminated by outer boundary effects. Refer- 
ring still to the state, we again have the situation where 
for a given time r and fixed value of H sp i n , the Im/i tends 
to have pervaded more of the superconductor region than 
the real part. The n amplitudes have some similarities 
with those in Fig. in that Re/i is smaller and has a 
weak dependence on the scattering strength. Although 
Im/i demonstrates a stronger dependence on H sp i n , and 
besides prominently peaking at the interfaces, the overall 
magnitudes for both components are effectively reduced 
for all values of r shown. 

This penetration of the triplet correlations into the S 
material can be conveniently described in terms of time- 
dependent penetration depths. These can be calculated 
for either /o or f\ . We will focus here on the real parts of 
/i, with the understanding that an analogous approach 
could be followed for either fa or the imaginary com- 
ponents. However, it is also evident from Fig. [2] that 
for fa, the penetration of triplet correlations into S is 
only weakly dependent on H sp i n , which is consistent with 
the fact that /o does not emerge from H sp i n only. The 
method to extract any sort of characteristic length de- 
pends of course on the problem at hand, and for the /i 
amplitudes we find the following definition yields sensible 
results, 



£(t) 



/ s dy|y-y ||Re{/i(r,r)}| 

f s dY\Be{fi(X,T)}\ 



(3.1) 



which is slightly modified from that used previously.— 
The definition Eq. (|3.ip accounts better for cases where 
the overall shape of the amplitudes (see Figs. [2] and [3]) 
varies depending on the parameter values. The coordi- 
nate shift, Y — Yq, accounts for measuring the distance 
from the interface: for our coordinates, Yq = Dp/2 and 
the integration extends over the S region. This defini- 
tion gives the expected result if the function f\ were a 
pure decaying exponential. The values of the dimension- 
less I are in units of kp. The results are plotted in 
Fig. |4] as a function of r. These results are for the same 
range of H sp i n as in the previous figures, as shown in 
the legend. Both the (top) and ir (bottom) state pen- 
etration depths are shown. There is an approximately 
linear behavior, in both and w cases, at earlier times 
and a deviation from linearity at later times. For the 
state, the values of I reach discernible maxima that 
get shifted to larger times with increased spin-dependent 
scattering rates. For larger values of H sp i n , the times nec- 
essary to reach the peak would presumably correspond to 
times where these triplet correlations would have reached 
the boundary and finite size effects would be a concern. 
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FIG. 4: (Color online) The penetration depth, as defined in 
dimensionless units in Eq. (|3.1[) . of the equal spin /i triplet 
component into the S material, plotted as a function of r for 
both the "0" and 'V configurations (as labeled). 



This increase for larger H sp i n , however, should not con- 
tinue much beyond the range shown, since we know-i£ 
that the triplet correlations do eventually begin to de- 
cay after a given characteristic time r of a few times 2ir. 
It is clear, however, that the penetration extends over a 
wide range of times over regions much larger than the 
superconducting coherence length. The results here are 
consistent with what we saw in Fig. [21 where at r = 20, 
we see for H sp i n = 0.2,0.4, and 0.6, the manifestation 
of secondary broad maxima at large Y, and little varia- 
tion among the f\ amplitudes, while the other values (0.8 
and 1) have generally a monotonic decline, and thus for 
these higher values they have smaller characteristic pen- 
etration depths. The penetration of equal spin triplet 
correlations into a tt junction is much less dependent on 
the spin scattering strength than for junctions except 
of course at very small H sp i n . This is again consistent 
with what was shown in Fig. [21 where the real parts of 
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the /i correlations demonstrated very little dependence 
on H spin , and were weaker overall, away from the inter- 
face, than their junction counterparts. This is reflected 
in the penetration depth behavior, where the depth is 
reduced compared to the 7r case and the penetration is 
similar for the broad range of spin dependent scattering 
strength. The only exception being the extreme case of 
H S pi n = 1, which always has a greater penetration, but 
this too agrees well with what is observed in Fig. [2] 

The density of states, measured in principle by in STM 
experiments, is one way of probing indirect evidence of 
the triplet superconducting states, and carries valuable 
energy resolved spectroscopic information. We therefore 
show next, in Fig. the DOS, computed from Eq. (|2.9[) . 
For computational purposes, we represent the delta func- 
tions in Eq. (|2.9[) as the low T limit of the derivative 
of the corresponding Fermi functions. We find that a 
fine mesh of e± is necessary to properly calculate the 
energy resolved DOS, and that there are significant con- 
tributions to the DOS from both longitudinal propaga- 
tion (small £jj and from large off- normal incidence at 
the interface (large e±). The junction here is therefore 
not appropriately described within the tunneling limit, 29 
and would yield differing results if only a narrow tun- 
neling cone was used in the calculation. By considering 
spin active interfaces, it was shown2£ that various signa- 
tures arise in the DOS, including a unique double gap 
structure. In Fig. [5] the DOS is shown, normalized to 
its bulk value in S, as a function of energy e (relative to 
Ep ) , in units of the bulk Ao . We consider four values of 
H sp in for both the state (left panels) and the it state 
(right panels). The top panels shows results summed 
over both spins and averaged over the entire thin F layer 
(label "F") while the bottom panels shows results (also 
summed over spins) averaged over the whole length of one 
of the superconductors (label "S"). For the S regions, in 
both the and n junctions there are BCS-type peaks 
at e/Ao « ±1, reflecting the bulk-like behavior. Inside 
the region of the bulk gap there is a secondary struc- 
ture reflecting Andreev states. These secondary peaks 
were also found in Ref. [3l|. We also see that the subgap 
peaks arise even in the absence of interface spin activity 
(at |e/Ao| w 0.5), and originate mainly from individual 
spin channels, depending on the sign of the energy: the 
prominent subgap peak at negative energy is due to the 
occupation of spin-up quasiparticles, N^, while its pos- 
itive energy counterpart arises from iVj_. The position 
of these subgap peaks varies with H sp i n in a way that 
seems to reflect the condensation energy in Fig. [TJ par- 
ticularly in the state. At higher energies, |e/Ao| > 1, 
both spin bands contribute equally to the DOS. In gen- 
eral within S and for the range of energies shown, the 
approximate relation, A|(e) ~ iVj(— e), holds, giving the 
observed symmetry in energy for the total DOS in the S 
region (Fig. [5] bottom panels). Increasing H spin , tends to 
flip the spins at the interface, and thus bound states in S 
predominantly occupied by a given spin species, become 
replaced by the opposite spin quasiparticles. This is con- 



firmed by examination of the individual spin density of 
states (Eq. 12. 9p . The half metallic ferromagnet modi- 
fies this bound state picture in that region due to the 
existence of only one spin band at the Fermi level (top 
panels). Here the majority of the energy-resolved states 
must come from spin-up quasiparticles. The spin-flip pro- 
cesses inherent to the scattering events at the interface 
(the parameters H sp i n or 7) in conjunction with proxim- 
ity effects can however cause an enhancement of subgap 
bound states at small energies attributed to a small num- 
ber of minority spin states in the magnet. Also due to 
the strong spin splitting in that region, there is signifi- 
cant particle-hole asymmetry. At higher values of |e/Ao| 
the correct limit is approached (~ (1/2)(1 + 1) 1 ^ 2 )- The 
structure in the region |e/Ao| < 1 is now considerably 
more complicated but consistent with what is seen in the 
S region: the number of states at zero energy is greater 
when the condensation energy is very small (see Fig. [I]). 

We turn now to the trends the peaks in S follow as a 
function of H S pin, for both junction configurations. As 
Hgpin is increased from to about 0.8, the value at which 
the condensation free energies are near their minimum 
values (see Fig.[lJ, we see that these peaks tend to merge 
and there is a zero-energy single or double peak signature, 
depending on whether it is a 7r or junction respectively. 
For Hgpin — 1) the peaks widen, to nearly the same en- 
ergies as for Hgpin =0.4. It does appear that for both 
junction states, depending on the spin strengths, there 
exists a small subgap region which resembles an energy 
gap in the DOS. Strictly speaking, however, the energy 
spectrum is gapless for the whole range from no spin flip- 
ping to strong activity at the interface; there is always a 
finite, albeit small in some cases, number of states within 
this Andreev bound state region. Adding a sufficiently 
strong spin-independent scattering or further increasing 
ds would eventually create a region with no states. The 
existence of the half metal and the interface scattering 
thus still influences the superconductor when consider- 
ing spatially averaged behavior. 

To gain further insight into the relative proximity ef- 
fects inherent to these junctions, we present in Fig. [5] 
results for the influence of H ap i n on the inverse proxim- 
ity effect, that is, on the local magnetic moment as de- 
fined in Eq. (|2.8p . We find that if the large e± off-normal 
trajectories are not fully included, the magnetic moment 
does not reach its proper normalized values and spatial 
properties. The calculation of m z and m x thus serves 
as another check to ensure that all the requisite states 
are included for other calculations. The results for the 
local magnetization (as defined above in Eq. (|2.8p and 
normalized by (1b) are shown in Fig. [6] They display the 
penetration of the magnetization into the superconductor 
region, as well as the weakening spin polarization in F. 
This is currently a topic of extreme interest, experimen- 
tally and theoretically, especially when trying to clarify 
the complicated spin structure in these systems. One can 
also view the introduction of spin scattering as an ultra 
narrow domain wall at the interfaces. In the figure both 
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FIG. 5: (Color online) The local density of states, normalized to unity in the normal state of the S material. The top panels 
show results within the thin F layer, while the bottom one is for the superconductor. The left column is for the state and 
the right one for the tv state. 



the x and the z components of the local magnetic moment 
are plotted as a function of Y for several values of H sp i n 
in the range < H spin < 1. The results are the same for 
and 7r states. Results for the state are shown. The 
x component vanishes by symmetry since the exchange 
field lies in the z direction when H sp i n — 0. At nonzero 
values of this parameter, m x grows very quickly in the 
interface region, while remaining zero in the center of the 
F layer and also, of course, deep in the superconductor. 
For this reason only the central part of the system is 
included in the plot. The z component (bottom panel) 
has a very weak variation with H spm near the center, 
Y = 0, of the half metal, but it does show a dependence 
on the scattering strength in the superconductor near 
the interface. The inverse proximity effect is clearly evi- 
dent where the induced magnetization component m z in 
the S region near the interface is oppositely directed to 
that in F, effectively screening the magnetization in F, by 



an amount that increases monotonically with H sp i n . On 
the ferromagnet side very near the interface, this compo- 
nent of the magnetic moment correspondingly weakens 
with increased H spini before rising up to near the half 
metallic bulk value of unity. The screening effect in the 
superconductor is apparently stronger for m z than m x . 
For the latter the induced magnetic moment is very sim- 
ilar in adjacent regions near the boundaries, and there 
is near symmetry about the interfaces. The component 
m x reverses signs in both F and S for all H sp i n , while 
m z is briefly negative for only the larger H sp i n , demon- 
strating the competing effects from the exchange field 
and spin scattering strength. The observed spatial char- 
acteristic of each component in Fig. [5] reveal that m in 
the vicinity of the interfaces tends to not only change 
magnitude as a function of position, but it also rotates. 
The magnetization also changes direction as a function 
of H sp i n for fixed Y, illustrating again the important role 
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the proximity effects play on the relevant self-consistent 
quasiparticle wavefunctions and energies. 



IV. CONCLUSIONS 

We have investigated the effect of interfacial spin-flip 
scattering on the triplet correlations that emerge in an 
SFS trilayer. We have studied this system by solving 
in a fully self-consistent way the BdG equations in the 
clean limit. We have considered both and it junctions 
and found that the results depend strongly on the junc- 
tion state. Triplet amplitudes, odd in time as required 
by the Pauli principle, have been found to exist and we 
have studied them in detail for the case where F is half- 
metallic. We have found that the m = ±1 triplet am- 
plitudes emerge and then subsequently increase (at finite 
times) very rapidly with the dimensionless spin flip pa- 
rameter Hspin ■ The degree at which the equal-spin triplet 
correlations pervade the S layers has been discussed in 
connection with the respective penetration depths. We 
also have presented results for the local energy resolved 
DOS averaged over both the S and F regions as a func- 
tion of the spin-flip rate. The or 7r state signatures may 
provide clues as to how different-time triplet states indi- 
rectly influence the subgap energy spectrum. We have 
also considered the inverse proximity effect (the penetra- 
tion of the magnetization into S and its weakening in F) 
and found that near the interfaces the magnetization ro- 
tates as a function of position or of H sp i n . Ultimately, the 
induced spin imbalance in the superconductor effectively 
screens the polarizing effects of the half metal. 



FIG. 6: (Color online) The x (top panel) and z components 
(bottom panel) of the local magnetic moment normalized by 
/is (see Eq. (12. 8[) ) plotted as a function of the dimensionless 
coordinate Y for H sp i n = from to 1 at 0.2 intervals. Only a 
limited range of Y is included. The and it state results are 
identical for the entire range of spin flipping considered here. 
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